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Abstract 



O ■ The 0+ states of 12 C are considered within the framework of the microscopic three-a-cluster 

^t- : 

model. The main attention is paid to accurate calculation of the width of the extremely narrow 



near-threshold 0^ state which plays a key role in stellar nucleosynthesis. It is shown that the 
0^-state decays by means of the sequential mechanism 12 C — > a + 8 Be — ► 3a. Calculations are 
performed for a number of effective a — a potentials which are chosen to reproduce both energy 
and width of 8 Be. The parameters of the additional three-body potential are chosen to fix both 
the ground and excited state energies at the experimental values. The dependence of the width 



5h 

on the parameters of the effective a — a potential is studied in order to impose restrictions on the 
potentials. 

PACS numbers: 21.45. +v, 21.60.Gx, 23.60.+e, 24.30.Gd 
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I. INTRODUCTION 



The processes with few (three and more) charged particles in the initial or final state 
are complicated phenomena which so far have not been completely understood. The main 
difficulty stems from the necessity to describe the continuum wave function of three (or 
more) charged particles (three-body continuum). Reliable description of the continuum 
three-body wave function is of importance for a number of problems in nuclear physics and 
nuclear astrophysics. As the first example one should mention the famous nuclear reaction 
- formation of the 12 C nucleus in the triple-a low-energy collisions. This reaction is of key 
importance for stellar nucleosynthesis as a unique possibility for helium burning that 

allows further synthesis of heavier elements. Other interesting examples of the three-body 
nuclear processes are double-proton radioactivity which has been a subject of thorough 
experimental and theoretical investigations during the last years (more details can be found 
in the recent reviews (3, Q|) and decay of the long-lived 1 + state of the 12 C nucleus j^J. 
For the problems of this kind, even qualitative understanding of the reaction mechanism 
is crucial. In this respect, Coulomb-correlated penetration of outgoing particles through a 
multidimensional potential barrier has been considered in Ref. ja| thus describing qualitative 
features of multicluster decay of atomic nuclei. 

Of key importance for description of the triple-a reaction are both the near-threshold 
three-body resonance (0^ state of 12 C) predicted in Ref. [J| as the only explanation for 
observable abundance of elements in the universe and the low-energy a — a resonance (the 
ground state of 8 Be). Due to existence of these resonances, sufficiently fast helium burning 
in stars was explained by the sequential mechanism 3a — > 8 Be + a — > 12 C(02~) — > 12 C + 7 
of the reaction. Indeed, the predicted 0jj~ state of the 12 C nucleus was observed in the 
experiments [3, 0| and was studied in the later works; in particular, the decay mechanism 
was a subject of investigation in Ref. jsj]. The experimental studies must be supplemented 
by microscopic calculations to provide unambiguous determination of the decay mechanism 
and resonance width V ~ 8eV the extremely small on the nuclear scale. 

Besides the resonance triple-a reaction, it is of interest to consider in astrophysical appli- 
cations, as pointed out in Ref. 10], the non-resonance reaction 3a — > 12 C which takes place 
at low temperatures and high densities. Helium burning at such conditions is possible in 
accretion on white dwarfs and neutron stars. The non-resonance reaction was considered in 
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a number of papers 111 12l, llii liii based on the model assumptions; however, a consistent 
treatment of the three-body dynamics is lacking. In this respect, note that at ultra- low en- 
ergies any approximation can lead to an error of a few orders of magnitude in the calculated 
reaction rate. 

Besides astrophysical applications, studies of the three-a scattering provides important 
information about the effective a — a interactions which is of interest for the a-cluster 
calculations. As the a-particle is the most tightly bound nucleus, many low-energy nuclear 
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properties can be successfully calculated within the framework of the a-cluster model 
fifiy Q. GeoeraH., the three-hody calculus allow one t o reduce the once rtamty 
in the two-body potential which can be hardly determined only from the two-body data. 
One of the principal opportunities for unambiguous determination of the a — a effective 
potential is to set the calculated width of the three-body resonance to its experimental 
value. In addition, one should mention that recently the near-threshold a-cluster states 
have attracted a special attention in connection with a-particle condensation in low-density 



nuclear matter 
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In this paper, properties of the + states of 12 C are considered using the 3a-cluster 
model with the main emphasis on the calculation of the width of the near-threshold 0^ 
state. Calculations are performed for a number of effective a — a potentials which are 
chosen to reproduce with a good accuracy both the energy and the width of the a — a 
resonance (ground state of 8 Be). Furthermore, due to the strong exponential dependence of 
the resonance width on the resonance energy, calculation of the width makes sense only if 
the resonance position is fixed. For the 3-a resonance under consideration, this requirement 
is satisfied by adjusting the parameters of the additional three-body potential which must 
be introduced to describe the effect of a-particle distortions near the triple-collision point. 
More precisely, the three-body potential is chosen to fix both the ground and excited state 
energies at the experimental values. Following Ref. [2^, the method of calculation is based 
on the expansion of the total wave function in terms of the eigenfunctions on a hypersphere 
(at a fixed hyper-radius), which allows solving both the eigenvalue problem for the ground 
state and the scattering problem for the excited resonance state. The eigenfunctions on 
a hypersphere are calculated by using the variational method with a flexible set of trial 
functions which describe properly the three-body wave function both at large and small 
interparticle distances. 
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The present three-body calculation of the near-threshold resonance is the first necessary 
step in the unified treatment of the low-energy triple-a reaction. Both the method and the 
numerical procedure can be used to calculate the reaction rate at lower energies where the 
resonance mechanism turns to the non-resonance one. 



II. METHOD 

The present paper is aimed at microscopic description of the low-energy scattering of 
three a-particles whose features are to much extent determined by the two- and three-body 
resonances. In this respect, the principal problem is reliable calculation of characteristics 
of the extremely narrow near-threshold resonance (O^-state of 12 C). The a-cluster model 
is used that allows for taking account of the most important features of the wave function, 
i.e., the three-a-cluster and two-cluster a + 8 Be components. All the effects connected with 
both the internal structure of a-particles and the identity of nucleons are incorporated in the 
effective a — a potentia 



form as in papers 
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_3esides, the additional three-body potential of a simple Gaussian 
171 ] is introduced to describe the effects beyond the three-cluster 
approximation. Thus, the model allows description of both the ground and the excited + 
states of 12 C. Considering the challenging problem of reliable calculation of the resonance 
width, the effective two-body potential should satisfy the restriction that the position and 
width of 8 Be are fixed at the experimental values. In a similar way, the three-body potential 
will be chosen to obtain experimental energies for both ground and excited states of 12 C. 

For the low-energy scattering in question, one should consider only the + states (the 
total angular momentum L = 0). The units h = m = e = 1 are used throughout the paper 
unless other is specified. The Schrodinger equation for three a-particles reads 

-A Xi - A y , + £ V(x 3 ) + V 3 (p) - E j = o (1) 

where the scaled Jacobi coordinates are Xj = — r k , yj = (2r^ — Tj — r k )/\/3 and is the 
position vector of the ith particle. In the following it is convenient to use the hyperspherical 
coordinates 0<p<oo,0<aj, and < $i < n defined as 

a i ■ a i n ( X »yQ (o \ 

Xi = pcos—, yi = psm—, cos^ = . (2) 

2 2 Xiyi 
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In the Schrodinger equation (JIJ the effective two-body potential V(x) is a sum of the short- 
range and Coulomb interactions 

A 

(3) 



V(x) = V s (x) + - 
x 



where the short-range a — a potential 



V„(x) = V r e-^ x2 - V a e~ Ax2 



(4) 



is obtained by modification of the Ali-Bodmer potentials [24]. The three-body potential 



V 3 (p) = V e'^ 2 



is chosen as the same function of the hyper-radius p as used in Refs (Ha, |1 



(5) 



A. Eigenfunctions on the hypersphere 



In terms of the hyperspherical variables the Schrodinger equation (P) for L = reads 
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(7) 



is the grand angular momentum operator up to a constant factor. In order to solve both 
the eigenvalue and scattering problems for Eq. (jSJ the total wave function is expanded in a 
series 

y = p- 5/2 J2fn(p)®n(n,6,p) (8) 

n 

on a discrete set of eigenfunctions $ n of the following equation on the hypersphere 

.2 3 



A*-^E^(pcos|)+A„(p) 



$ n (a,0,p) = , 



(9) 



as proposed in Ref. [2^|. At each p the index n = 1,2,3, .. . enumerates the eigenvalues 
A n in ascending order and the eigenfunctions <& n (a,9,p) are normalized by the conditions 
(<& n |<& m ) = S nm where the notation (-|-) means the integration over the invariant volume 
on the hypersphere dfl = sin 2 aidctid cos 6^. Due to identity of a-particles both the total 
wave function \1/ and the eigenfunctions Q n (a,8,p) are symmetric under any permutation 
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of particles i,j, and k. Given the expansion (JHJ) of the total wave function, the Schrodinger 
equation © is reduced to the system of hyper-radial equations (HRE) 



Up) + 



dp 2 p 2 \ 4 
( d d \ 

Yl \ Qnm(p)g^ + -Q- p Qnm{p) ~ P nm (p) fm(p) = (10) 



where 



Qmn(p) = (®rn 

/d$ m 



Pmn (p) 



dp 



d$ n \ 
dp J 

dp 



(11) 
(12) 



The coefficients X n (p), Q nm (p), and P nm (p) of HRE (fTU|) are calculated using the vari- 
ational method for solution of the eigenvalue problem Q. The variational basis consists 
of TV trial functions Xi with the same symmetry under permutations of particles as the 
eigenfunctions Q n (a,8, p). In view of an essentially different structure of the eigenf unctions 
$ n (a,#, p) at different values of p, it is necessary to use a flexible basis of trial functions 
which allows one to describe the two- and three-cluster structure of the wave function in the 
asymptotic region. 

First of all, the basis contains a set of the symmetric hyperspherical harmonics (SHH) 
H nm which are the eigenfunctions of the operator A*, i. e., 

[A* + K(K + 2)] H nm = (13) 

where K = In + 3m, the non-negative numbers n and m enumerate SHH, and 2K is the 
order of SHH. For explicit construction of SHH it is convenient to use another set of the 
hyperspherical variables < £ < — , — 7r < <Pi < 7r j^Ql defined by 

sin£ = sin a j sin ^j, 

cos£ cos ifi = cosctj, (14) 
cos £ sin ifi = sin cos 9i . 

In these variables 

H n mtt,v) ~cos 3m ^Pi ' 3m )(cos20T 3m (cos^) ~ 3™ (2£) cos 3m<p (15) 

2 ' 2 
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where P^ a '^(x) and T n (x) are the Jacobi and Chebyshev polynomials and d 3 mk (j3) is the 
Wigner function. The variable £ is invariant under permutations of particles and, therefore, 
is independent of the index i enumerating the Jacobi variables. On the other hand, (pi 
changes to p>i ± 2tt/3 under the cyclic permutations as | tpi — ipj \= 2tt/3 and ifi — > —<fi 
under the permutation of particles j and k. As follows from Eq. (fT3j) and the above properties 
of the variables £ and ipi, SHH are completely symmetric under any permutation. 

In the numerical calculations, the basis of trial functions contains a set of all SHH 
0) = if njmj (£, tp) with those indices rij and for which K does not exceed the maxi- 
mum value K max , i. e., iQ = 2rii + 3mj < K max . One can count that the total number of 
such SHH for which 2rij + 3mj < K equals K(K + 6)/12 + 1 for K being a multiple of 6 and 
([if/6] + 1)(K — 3[if/6]) otherwise. Here [x] stands for the entire part of x. Usage of SHH 
in the basis of trial functions provides an excellent description of the eigenf unctions at small 
p, where the kinetic energy term dominates, and quite a good description at intermediate p, 
where the cluster effects still do not dominate. However, the two-cluster component of the 
wave function corresponding to the configuration a + 8 Be can be hardly described by a set 
of SHH due to rather slow convergence that hinders the calculation at sufficiently large p. 

In order to describe the two-cluster configuration, the basis of trial functions should also 
include the p-dependent symmetric combinations 

3 a- 

Xi {a, B) = Yj <Mp cos ( 16 ) 
i=i 

of the two-body functions <fii(x) which are chosen to describe the wave function of the two- 
body a — a resonance. More precisely, a set of <pi(x) includes a few Gaussian functions 

4>i(x) = exp (-fax 2 ) (17) 

which allows the two-body wave function to be described, with properly chosen parameters 
within the range of the nuclear potential V s (r). In addition, the function 

<f,{x) = x 1/4 exp (-4y^(l + ax)) (18) 

is used to describe the two-body wave function in the under-barrier region. This latter 
function is of the asymptotic form of the Coulomb wave function which is optionally cut off 
by the parameter a at large distances. 
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Although the eigenvalues A n (p) are directly determined in the variational calculation, the 

coupling terms Q nm (p) and P nm {p) can be hardly determined by means of definitions (fTTj) 

and (fP2l . which is hindered due to necessity to calculate the derivatives For this 

op 

reason, Q nm (p) are calculated by using the exact expression 



3 . x _! / T q , a 2 <9Vs(pcosf) 



Qmn(p) = | (An ~ A m ) 1 ^$ 



+ 2pK(pcos-) +p 2 



$n ) (19) 



cos a 2 dp 

which is derived by differentiating the eigenvalue equation (J7J) with respect to p and project- 
ing the result on the function $ m . Furthermore, P nm (p) are calculated by using the exact 
sum rule P = — Q 2 for the matrices P and Q, which leads to the approximation 

N 

Pmn(p) = Y.Qml(p)Qnl(p) (20) 
1=1 

on the limited basis of N trial functions. 



B. Boundary conditions and characteristics of C states 

Properties of the ground 0i~ state and the excited 0^ resonance are determined by solving 
the eigenvalue problem (at E < 0) and scattering problem (at E > 0) for HRE ()10j). 
respectively. Denote the hyper-radial functions as f^\p) f° r the ground state and f^\p) 
for the scattering problem at energy E. According to (jHJ), all these functions satisfy the zero 
boundary conditions at p = 0. The square integrable solution of HRE (fTUjl satisfying the 
condition 

oo 

I \ti\p)\ 2 dp = 1 (21) 
n J Q 

unambiguously determines the energy E gs and the wave function of the ground 0i~ state. 

The position E r and width T of the near-threshold 0^ resonance are calculated by solving 

HRE ()10|) with the asymptotic boundary conditions corresponding to the ingoing wave in 

the first channel, i. e., ff(p) is a sum of the ingoing and outgoing waves in the effective 

potential U\(p) = — [4Ai(p) + 15/4]. More precisely, the asymptotic boundary conditions 
P 

are imposed near the turning point p t of the first-channel effective potential defined by the 
condition Ui(p t ) = E. As expected, calculations reveal (see also Fig. [TJ) that at p ~ p t the 
effective potential Ui(p) is to a good approximation expressed as 

U 1 (p)^E 2a + ^ (22) 



where Em is the energy of the two-body resonance (the ground state of 8 Be) and the Coulomb 
parameter q = 16/v3- In fact, r.h.s. of Eq. (}2*2"|) is the energy of the two-cluster system 
a + 8 Be at a large fixed hyper-radius p. For the scattering at energy above the two-body 
resonance (E > i?2«), in view of expression (f2*2*|) . the first channel hyper-radial function can 
be written as 

f[ E \p) ~ F ( V , kp) + tmS E G { V , kp) (23) 
in the range of hyper-radius values p ~ p t . Here the wave number in the first channel 



k = y/E — E2a, Fo(r],kp) and Go(rj,kp) are the Coulomb functions with the parameter 
i] = 8/(v / 3^), and 5e is the scattering phase shift. Due to strong repulsive potentials 
U n {p) for n > 2 the outgoing waves in the upper channels are negligible at small energies 
E < lMeV. This allows the zero boundary conditions 

/f \P) = (24) 

to be imposed at some value of the hyper-radius p > pt for all n > 2. The resonance 
position E r and width T as well as the non-resonant phase shift Sb g are defined by fitting 
the calculated near- resonance phase shift 5e to the Wigner dependence on energy 

cot(5 E - 5 bg ) = ^(E - E r ) . (25) 

In the following description of the 0^ state it is suitable to treat the ultra-narrow res- 
onance as a true bound state with the wave function f^'(p) ~ fi Er Kp) corresponding to 
the scattering solution at the resonance energy E r and normalized on the finite interval 
< p < Pt by the condition 

Y,J\fi 2) (p)\ 2 dp=l. (26) 

n J Q 

It is of interest to determine also the root-mean-square (RMS) radii 

i N t 

R{1) = W E(* W ( r * " R -) 2 (27) 
iV < k 

of the ground (i — 1) and excited (i = 2) states and the monopole transition matrix element 

N p 

M 12 = ^«|(r fc -R cm ) 2 |^) . (28) 

k 
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A sum is taken over N t nucleons in (|27J1 and over N p protons in (|28|) and R cm is the center-of- 
mass position vector. Following the definitions (J2*7|) and (|2*Sjl . in the three-a-particle model 
one obtains the expressions 

r {i) = y^+|S . (29) 

Pi 

M 12 = Y.[fL 2) (p)rt\p)P 2 dP (30) 

where i? a = 1.47fm is the RMS radius of the a-particle and the RMS value of the hyper- 
radius pi for the ith. state in the three-body model is defined by 



Pi 



Y.J\ff(p)\ 2 p 2 dp. (31) 



III. NUMERICAL RESULTS 



Calculations have been performed with four a — a potentials which are obtained by 



modification of potentials (a) and (d) from Ref. |24j. The parameters of the potentials 
have been chosen to reproduce the experimental value of the a — a resonance ( 8 Be) energy 
E>2u = 91.89keV whereas resonance widths have been allowed to vary within the experimental 
uncertainty ±1.7eV. The parameters of the potentials and the corresponding a — a resonance 
widths 7 are presented in Table HJ Only the strengths V r and V a of the repulsive and attractive 
parts are varied for potentials 1-3, while the parameters \L r and \i a are modified for "harder" 
potential 4. This choice of the potentials makes it possible to study the dependence of the 

TABLE I: Parameters of the short-range a — a potential © and the corresponding widths 7 of 
the a — a resonance (the ground state of 8 Be). 

V r (MeV) Mfm- 1 ) V a (MeV) fi a (fin" 1 ) 7 (eV) 



1 


82.563 


1/1.53 


26.1 


1/2.85 


6.80 


2 


279.206 


1/1.53 


40 


1/2.85 


8.53 


3 


20.012 


1/1.53 


16.5 


1/2.85 


5.11 


4 


197.680 


0.7 


80 


0.475 


5.10 



three-body characteristics on the shape of the two-body potential. 

The energies of the ground and excited states E gs and E r , width of the excited state T, 
RMS radii and monopole transition matrix element M 12 for the three-a system are 
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calculated, as discussed in the preceding section, by numerical solution of HRE (fTUj) with 
boundary conditions (}2*3*)) and (|24|) . For the eigenvalue problem, numerical integration in the 
interval < p < 25fm provides the relative accuracy not worse than 10 -6 for the calculated 
binding energy. For the scattering problem, numerical integration is carried out in the 
interval < p < p max so that p max ~ 50fm is chosen beyond the turning point p t ~ 45fm 
for the first-channel effective potential U\{p). Choice of the parameters of the numerical 
integration does not affect the final results, i. e., the accuracy of the calculated E gs , E r , T, 
F&\ M 12 depends merely on the accuracy of the numerical calculation for A n (p), P mn {p), 
and Qmn{p)- In every calculation, i. e., for each two-body potential and each number of 
HRE Ni, the parameters Vq and b of the three-body potential V%(p) are chosen to fix the 
calculated E r and E gs at the experimental values E r = 0.3795MeV and E gs = — 7.2747MeV. 

For the calculation of the HRE coefficients A„(p), P m n{p), and Q mn {p) by the variational 
method, two different sets of trial functions were used in the interval p < 20fm and in 
the asymptotic region p > 20fm. At p < 20fm the basis of trial functions contains 147 
SHH, which corresponds to K max = 39. Recall that convergence with a number of SHH 
becomes very slow with increasing hyper-radius. For this reason, at p > 20fm the basis of 
trial functions contains 108 SHH (K max = 33) and four trial functions which describe the 
cluster configuration a + 8 Be, namely, three functions of the form (fTTj) and one function 
of the form ()18|). The p-independent parameters $ in functions (fTTj) were determined by 
minimizing the first eigenvalue Ai(p). The calculated A n (p), P m „(p), and Q mn (p) practically 
do not depend on the basis of trial functions near p = 20fm, which allows matching the 
results of calculation with two different basises. The three lowest effective potentials U n (p) = 
— [4A„(p) + 15/4] + P nn (p) (n = 1 — 3) for two-body potential 3 are depicted for illustration 
in Fig. [TJ 

The detailed information on the three-body system can be obtained by considering the 
eigenfunctions on a hypersphere <p, p). Due to symmetry, $„(£,<p, p) are periodic 

functions of the variable <p with the period 27r/3. As the main contribution to the total wave 
function comes from the first term (n = 1) of expansion (jHJ), the first eigenfunction $i(£, <p, p) 
practically determines the structure of the system. For illustration, three-dimensional plots 
of the first eigenfunction <3>!(£,(p,p) (for two-body potential 3) are shown in Fig. [2 at three 
values of the hyper-radius p = 5fm, p = 15fm, and p = 45fm that correspond, as seen in 
Fig. to the minimum, the maximum, and the turning point of Ui(p) at E = 0.3795MeV. 



11 



FIG. 1: Three lowest effective potentials U n (p). Shown are also the asymptotic two-cluster depen- 
dence U\(p) « E2 a + - (thin dashed line) and the energy of the 0^" state E = 0.3795MeV (thin 
horizontal line). 

U n 
4 

2 



-2 

-4 

20 40 

P 

For all values of the hyper-radius, $i(£, (p, p) is small at the zero distance between a pair 
of a-particles, i. e., at the point £ = 0, <p = n/2 at the hypersphere, because of the strong 
repulsive term in the short-range potential V s . At sufficiently large p (near the turning 
point), $i(£, <p, p) exhibits a specific structure concentrated in the region where the two-body 
potential is attractive, i. e., around the point £ = 0, tp = tt/2. This structure corresponds to 
the two-cluster configuration a + 8 Be. For smaller values of p (as seen in Fig. [21 at p = 15fm) 
the two-cluster structure widens; besides, visible values of $i(£, <p, p) appear both at the 
point £ = tt/2, which corresponds to the configuration of the equilateral triangle, and near 
the point £ = 0, p = (or £ = 0, <p = 2n/3), which corresponds to the linear configuration. 
Next, at small p near the minimum of Ui(p), the most important is the triangle configuration 
with a noticeable weight of the linear configuration and without any trace of the two-cluster 
structure. 

The accuracy of the calculation is estimated by observing the convergence with increasing 
number of SHH. As a result, the accuracy of the most significant effective potential U\{p) 
turns out to be not worse than leV in the entire interval < p < 20fm and is much 
better for smaller p. In a similar way, the relative accuracy of P m n{p) and Q m n(p) is better 
than 10~ 5 in the same interval < p < 20fm. As the hyper-radius increases beyond 
p = 20fm, the accuracy of the calculation decreases due to a more complicated structure of 
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FIG. 2: The first eigenfunction $1 p) at different p. 




the eigenfunctions $i(£, <p, p). Nevertheless, as shown in Fig. the effective potential Ui(p) 
is in good agreement with the asymptotic dependence (|22jl in the interval 40 < p < 60, 
thus pointing to the sufficiently high accuracy. In particular, for all the potentials used, a 
fit of Ui(p) to Eq. (|22|) gives the values of the Coulomb parameter q which differ from the 
expected value q = 13.3MeV by less than 0.14MeV-fm. The fitted values of Ei a differ from 
the experimental energy of 8 Be by less than 0.004MeV for potentials 1-3 and by about 
0.009MeV for potential 4. 

Note that representation (|2*2*|) of the effective potential Ui(p) in terms of the energy of 
the two-cluster system a + 8 Be confirms the sequential mechanism of O2" state decay with 
formation of a + 8 Be at the first step. Also, this conclusion follows from the genuine two- 
cluster form of the first-channel eigenfunction p) at p ~ pt- As shown in Fig. E]at 
p = 45fm, $i(£;0>p) practically coincides with the symmetric combination of the 8 Be wave 
functions. Thus the total width is determined by the two-cluster decay width. 

The calculated T, (i = 1,2), M 12 , and the parameters of the three-body potential Vq 
and b are presented in Table |H] for four two-body potentials and different number of HRE 
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Ni. These results are compared with the experimental values and calculations 

HQ (the 



TABLE II: Characteristics of the + states of three a-particles for four a-a potentials calculated 
with N% HRE. The width of the excited state T, the RMS radii and R^ 2 \ the monopole 



transition matrix element M12, anc 
The results of the calculations 16 



the parameters Vq and b of the three-body potential are given. 



17j | and the experimental values are given in the last rows. 





Nx 


r(eV) 


^(^(fm) 


i?( 2 )(fm) 


M12 (fm 2 


) Vb(MeV) 


6(fm) 




1 


25 


2.56 


4.1 


9.52 


-21.998 


4.8993 


1 


2 


20 


2.55 


4.1 


9.08 


-23.993 


4.6608 




3 


19 


2.55 


4.0 


9.03 


-24.048 


4.6530 




1 


37 


2.82 


4.4 


10.3 


-29.291 


5.1213 


2 


2 


27 


2.77 


4.2 


9.00 


-39.124 


4.5223 




3 


26 


2.77 


4.2 


8.85 


-39.947 


4.4858 




1 


13 


2.25 


3.6 


8.23 


-15.812 


4.4506 


3 


2 


11 


2.24 


3.6 


8.15 


-16.059 


4.3964 




3 


10 


2.24 


3.5 


8.14 


-16.066 


4.3942 




1 


21 


2.51 


4.1 


8.74 


-14.548 


5.7248 


4 


2 


16 


2.49 


4.0 


8.36 


-15.129 


5.4595 




3 


15 


2.49 


4.0 


8.21 


-15.285 


5.4118 


[16] 




20 


2.36 




6.54 


-96.8 


3.9/V2 


[12 




1300 


2.47 




8.36 


-23.32 


3.795-^/2 


Exp. 




8.5±1.0 


2.47 




5.7 







factor \pl in the last column appears due to different definitions of p in these papers). As 
shown in Table |Hj the convergence in a number of HRE is sufficiently fast and solution of 
three HRE (N\ = 3) allows the resonance width to be determined with the accuracy about 
leV. The final accuracy of T depends mainly on the accuracy of the variational calculation 
and can be estimated as a few eV. The values of the parameter b presented in Table |H] are 
quite reasonable since they are in agreement with the value p = 2\^2R a w 4.16fm, which 
corresponds to triple collision of three hard spheres with radii R a . 
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IV. DISCUSSION AND CONCLUSION 



The main result of the present calculation is accurate determination of an extremely 
narrow width T of the three-body resonance (the excited state of 12 C). As shown in 
Tables E and |HJ the width T of the three-body resonance depends on the underlying two- 
body a — a potential and essentially increases with increasing width 7 of the two-body 
a — a resonance. In fact, variation of 7 within the limits of the experimental uncertainty 
gives rise to a change of T by a factor 2.6; therefore, the precise value of 7 is necessary for 
determination of T. Note that for all the potentials considered the calculated T overestimate 
the experimental value T exp = 8.5 ± leV by a factor 1.2 — 3.1. However, the difference of 
the calculated width for potential 3 and the experimental value is of order of the theoretical 
and experimental uncertainties. In addition, comparing the results for potentials 3 and 
4, one may conclude that width T of the three-body resonance depends on the potential 
shape (on the parameters p r and /i a ). It should be emphasized that dependence of T on the 
parameters of the a — a interaction is rather complicated due to addition of the three-body 
potential V^{p) which is chosen to fix E gs and E r at the experimental values. Nevertheless, 
the effect of V^{p) is not overwhelming since V is mainly determined by penetrability of the 
potential barrier and all the three-body potentials used in the present calculation are rapidly 
decreasing with increasing p in the barrier region 10fm< p < 45fm. 

Similar dependence of 7 is revealed for the ground-state RMS radius It is interesting 
that the calculated = 2.55fm is close to the experimental value R^J P = 2.47fm for 
potential 1, for which the two-body width 7 coincides with the most probable experimental 
value j exp = 6.8eV. Similar to T, the calculated R^ essentially depends on the parameters 
p r and p a - Both the RMS radius R^ of the exited states and the monopole transition 
matrix element M 12 are weakly dependent on the parameters of the two-body potential. For 
all the potentials, the calculated values of M 12 significantly overestimate the experimental 
value 5.7fm 2 that clearly deserves further investigation. 

Comparison with the previous microscopic calculation shows that the accuracy of the 
present calculation is much better than in Ref. r| whereas the methods of calculation are 
similar to each other. Note that the a — a potential in Ref. jl^] was chosen to fix the RMS 
radius of the ground state at the experimental value R^ p = 2.47fm. With this potential, 
the energy of the three-body resonance is misplaced by 0.47MeV. This is essentially above 
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the experimental values and leads to an unreliably large resonance width. 

The above discussion allows the conclusion that calculation of three-body observables can 
be used to impose restriction on the effective two-body a — a potential. In the future one 
should look for a possibility of reproducing the experimental values by a refined choice of 
potentials, in particular, by using more complicated three-body potentials. 

In conclusion, the three-a-cluster model is used to calculate the characteristics of the + 
states in the 12 C nucleus. In particular, the width of the extremely narrow threshold 
state is calculated with a good accuracy. The dependence of the width on the parameters 
of the a — a potential is studied. It is proposed to use the calculation of the width for 
selection of a — a potentials. It is directly shown in the three-body calculation that the 
state decays by means of the sequential mechanism ( 12 C — > a + 8 Be — > 3a). This conclusion 
is in agreement with the experiment in which the rate of the direct decay 12 C — > 3a 
is estimated to be less than 1% of the total rate. The calculation of the near-threshold 
3a-resonance can be considered as a part of the general study of both resonant and non- 
resonant reaction 3a — > 12 C at low energy. The present approach is promising for calculation 
of the triple-a reaction rate at low energy below the three-body resonance. This provides 
an opportunity for unified treatment of the crossover from the resonant to the non-resonant 
mechanism of the reaction. 
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